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Abstract 

This  paper  is  motivated  by  the  diverse  array  of  border  threats,  ranging  from  terrorists  to 
arms  dealers  and  human  traffickers.  We  consider  a  moving  sensor  that  patrols  a  certain 
section  of  a  border  with  the  objective  to  detect  infiltrators  who  attempt  to  penetrate  that 
section.  Infiltrators  arrive  according  to  a  Poisson  process  along  the  border  with  a  specified 
distribution  of  arrival  location,  and  disappear  a  random  amount  of  time  after  their  arrival. 
The  measures  of  effectiveness  are  the  target  (infiltrator)  detection  rate  and  the  time  elapsed 
from  target  arrival  to  target  detection.  We  study  two  types  of  sensor  trajectories  that  have 
constant  endpoints,  are  periodic,  and  maintain  constant  speed:  (1)  a  sensor  that  jumps 
instantaneously  from  the  endpoint  back  to  the  starting-point,  and  (2)  a  sensor  that  moves 
continuously  back  and  forth.  The  controlled  parameters  (decision  variables)  are  the  starting 
and  end  points  of  the  patrolled  sector  and  the  velocity  of  the  sensor.  General  properties  of 
these  trajectories  are  investigated. 


1  Introduction 

The  problem  of  illegal  cross-border  activity,  in  particular  along  the  U.S. -Mexico  border,  has 
drawn  much  attention  in  recent  years  [13],  [7].  In  addition  to  people  who  seek  jobs  and  better 
life  in  the  U.S.,  illegal  border  crossing  has  also  been  associated  with  criminal  activities  such 
as  smuggling  and  drug  trafficking.  During  the  year  2000,  over  1  million  pounds  of  marijuana, 
30,000  pounds  of  cocaine  and  200  pounds  of  heroin  were  seized  along  the  U.S.-Mexico  border 
[6] .  Thus,  illegal  border-crossing  has  been  a  cause  for  much  concern  to  local,  state  and  federal 
law  enforcement  authorities.  Moreover,  while  so  far  there  has  not  been  any  definite  evidence 
of  terrorists  involvement  associated  with  this  illegal  activity,  the  concern  is  that  terrorists 
with  weapons  of  mass  destruction  could  possibly  penetrate  the  porous  border  in  much  the 
same  way  that  illegal  migrants  and  smugglers  do.  This  scenario  poses  a  serious  national 
security  threat  that  must  be  handled  by  the  authorities.  While  several  stakeholders  at  all 
levels  of  government  take  part  in  securing  the  U.S.  against  such  possible  events,  the  most 
important  stakeholder  in  border  security  is  the  U.S.  Border  Patrol,  whose  current  primary 
mission  is  to  prevent  terrorists  and  terrorist  weapons  from  entering  the  country  [2].  This 
paper  is  motivated  by  that  mission. 
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Concerning  border  security,  an  increasingly  common  approach  to  border  surveillance  is 
to  use  a  combination  of  long  range  and  short  range  sensors  [1] .  Generally,  there  is  a  tradeoff 
between  the  range  of  the  sensor  and  its  resolution;  long  range  sensors  provide  less  resolution 
than  short  range  sensors  [1].  Although  the  models  developed  in  this  paper  can  be  applied 
to  a  variety  of  surveillance  systems,  we  focus  on  two  systems:  a  mounted  scanning  camera 
that  scans  a  certain  sector  of  the  border  and  then  leaps  back  to  the  starting  point,  and  a 
UAV- mounted  sensor  that  flies  back  and  forth  between  the  two  endpoints  of  the  sector. 

There  are  two  different  types  of  UAVs:  drones  and  remotely  piloted  vehicles  (RPVs). 
Both  are  pilot-less,  but  drones  are  programmed  for  autonomous  flight,  whereas  RPVs  are 
actively  flown,  remotely,  by  a  ground  control  operator.  UAVs  have  played  key  roles  in  recent 
conflicts  [1],  providing  reconnaissance,  surveillance,  target  acquisition,  search  and  rescue, 
battle  damage  assessments  and  attack  capabilities.  An  example  of  such  a  system  that  is 
increasingly  present  in  U.S.  border  surveillance  activities  is  the  Predator  B  UAV  system, 
which  can  provide  long-endurance  surveillance  and  communications  relay  [5].  While  UAVs 
are  used  to  extend  the  surveillance  range,  especially  in  the  absence  of  ground  infrastructure, 
ground-based  electro-optical  sensors  achieve  greater  resolution.  A  variety  of  electro-optical 
imaging  sensors  are  employed  in  border  surveillance.  A  few  examples  of  the  technologies 
in  use  are  high-resolution  imaging,  motion  detection,  temperature-differentiation  and  night- 
vision  devices.  In  addition  there  is  wide  use  of  commercially  available  cameras  adapted  for 
daytime  surveillance,  infrared  video  (IR)  detection  systems  and  laser  illumination  systems 
that  enable  high-resolution  imaging.  Other  optical  components  include  computer-operated 
pan/tilt/zoom  cameras,  visible  or  near-infrared  illuminators  for  night  vision  with  conven¬ 
tional  cameras  and  image-intensihers  for  long-range  night  vision  [5].  The  combination  of 
different  sensors  mitigates  the  shortcomings  of  each  single  sensor.  For  example,  IR  provides 
detection  in  atmospheric  conditions  where  video  imaging  is  ineffective,  even  though  the  reso¬ 
lution  it  provides  is  relatively  low.  Laser  illumination  sensors  provide  a  combination  of  long 
range  and  high  resolution  images  [5]. 

Previous  work  on  linear  patrols  mainly  concerns  anti-submarine  warfare.  The  classic 
WWII  book  by  Koopman  [9]  lays  the  foundations  for  analyzing  the  performance  of  barrier 
patrols.  Washburn,  [18]  and  [19],  extends  the  basic  barrier  patrol  model.  Models  of  barrier 
patrol  have  focused  specihcally  on  anti-submarine  warfare,  in  which  the  acoustic  degradation 
of  the  sonar  sensor,  caused  by  increasing  the  patroller’s  speed,  is  incorporated  into  the 
model,  and  the  optimal  patrol  speed  is  found.  A  game-theoretic  analysis  of  choice  of  speeds 
by  inhltrator  and  patroller  originated  in  [8]  and  [12],  and  was  developed  analytically  for  a 
continuum  of  speed  choices  in  [10].  While  these  are  important  references,  we  are  not  aware 
of  any  past  or  current  research  that  is  applicable  to  our  problem  of  perimeter  protection. 
Although  no  queueing  occurs  in  our  models,  some  of  the  techniques  we  use  resemble  those 
employed  in  models  based  on  spatial  queues;  see,  for  example,  [3]. 

This  paper  sets  forth  models  for  a  moving  sensor  that  patrols  a  certain  sector  of  a  border. 
The  goal  is  to  hnd  a  patrol  policy  that  minimizes  the  inhltration  rate  across  that  sector 
(maximize  detections)  when  the  sensor  follows  a  periodic  trajectory  at  constant  speed.  The 
controlled  parameters  (decision  variables)  of  the  trajectory  are  the  beginning  and  end  points 
of  the  patrolled  sector  and  the  sensor’s  velocity.  We  view  the  main  contributions  of  this 
paper  as  follows.  First,  we  develop  models  of  border  surveillance  that  capture  the  situation 
where  inhltrators  may  disappear  before  sensor  detection,  and  therefore  traverse  the  border 
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undetected.  Second,  we  analyze  and  prove  several  properties  of  the  optimal  trajectories  that 
are  operationally  signihcant.  Finally,  we  provide  a  simulation  analysis  for  the  case  where 
identifying  the  inhltrators  is  not  instantaneous. 

The  type  of  surveillance  systems  described  in  this  paper  have  recently  been  brought  to 
public  attention  (see,  for  example,  [17])  by  the  Department  of  Homeland  Security,  which 
signed  a  contract  with  the  Boeing  Corp.  to  secure  the  U.S.  southern  border.  As  part  of  this 
contract,  tower-mounted  and  UAV-mounted  sensors  will  be  installed  on  the  U.S.-Mexico 
border  along  parts  of  the  state  of  Arizona. 

The  rest  of  the  paper  is  organized  as  follows.  Section  2  presents  the  basic  framework  of 
our  models.  Section  3  discusses  properties  of  the  leap-to-origin  trajectory,  and  section  4  deals 
with  the  back-and-forth  trajectory  (see  Figure  1).  Section  5  extends  the  baseline  scenario 
by  considering  the  effect  of  sensor  velocity  and  the  case  where  it  takes  time  to  investigate 
each  detected  target. 

2  Model  Preliminaries 

The  two  main  components  of  the  problem  are  a  sensor  and  targets  (inhltrators).  We  initially 
assume  that  the  sensor  is  perfect  and  the  detection  and  identihcation  process  is  instantaneous. 
In  particular,  the  detection  probability  is  independent  of  the  sensor  velocity.  We  set  the 
velocity,  without  loss  of  generality,  to  be  1.  These  assumptions  are  relaxed  later  on.  The 
sensor  scans  a  certain  sector  of  the  border  continuously  and  repetitively  and  our  objective 
is  to  determine  the  endpoints  of  that  sector — its  origin  and  destination  points.  We  assume 
that  these  points  do  not  change  during  the  surveillance  mission.  Two  consecutive  visits  to 
the  origin  determine  a  cycle  in  the  surveillance  process.  Targets  arrive  to  the  perimeter 
according  to  a  homogeneous  Poisson  process  with  rate  a.  More  specihcally,  let  {Ti  :  i  >  1) 
be  the  collection  of  arrival  epochs,  then  the  number  of  arrivals  by  the  time  t,t  >  0  is  given 
by 

CXD 

i=l 

where  /(■)  is  the  indicator  function. 

The  arrival  locations  of  the  targets  are  described  by  a  collection  of  independent  and  iden¬ 
tically  distributed  (i.i.d.)  continuous  random  variables  (X*  :  i  >  1)  with  common  probability 
distribution  Fx-  After  spending  a  certain  time  at  the  border  area,  a  target  may  successfully 
inhltrate  the  border  and  therefore  will  no  longer  be  available  for  detection.  The  period  of 
time  elapsed  from  target  arrival  until  inhltration  is  the  reneging  time  and  is  described  by 
a  collection  [Ri  :  i  >  1)  of  i.i.d.  non-negative  continuous  random  variables  with  common 
distribution  function  Fr.  Because  targets  may  renege,  the  problem  is  to  hnd  the  endpoints 
of  the  surveillance  sector  that  represent  the  best  tradeoff  between  two  types  of  undetected 
(lost)  targets:  those  who  are  outside  the  surveillance  sector  and  those  who  renege  within  it. 

The  following  assumption  is  satished  in  all  practical  situations  of  interest,  and  greatly 
simplihes  the  mathematical  treatment. 

Assumption  1  Throughout  the  paper  we  assume  that  random  variables  (Xj  :  i  >  1)  posses 
a  density  fx  that  is  unimodal  and  continuously  differentiable  on  its  support.  We  also  assume 
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Figure  1;  The  LO  trajectory  (left)  and  the  BF  trajectory  (right):  Straight  lines  represent 
sensor’s  trajectories;  the  signs  “x”,  and  “o”  represent  arrival,  detection  and  reneging  of 
a  target,  respectively.  Cycle  length  is  u  in  the  left  hgure,  and  2u  in  the  right  hgure. 

that  Fji  is  continuously  differentiable  on  its  support. 

We  consider  two  types  of  generic  sensors:  A  human-controlled  remote  video  camera  and 
a  UAV-mounted  sensor.  A  human-controlled  remote  camera  scans  the  sector  of  the  border 
[5]  and  is  typically  installed  at  the  top  of  a  tall  tower.  The  camera  starts  the  surveillance 
process  at  one  end  of  the  sector,  moves  at  constant  speed  to  a  certain  destination  point, 
and  then  leaps  back  to  the  origin  instantaneously.  This  is  referred  to  as  a  leap -to -origin 
(LO)  trajectory.  A  UAV-mounted  sensor  travels  along  the  border  sector  back  and  forth  in  a 
constant  speed,  in  a  back- and- forth  trajectory  (BF).  Figure  1  illustrates  the  LO  scenario  on 
the  left-hand  side  and  the  BF  scenario  on  the  right-hand  side. 

The  main  component  in  our  formulation  is  the  process  {Nci{t)  :  f  >  0),  which  counts  the 
number  of  detections  by  time  t.  Our  construction  guarantees  that  is  a  Poisson  process 
for  any  cycle  length  u,  with  rate  that  depends  on  Fx  and  Fr,  on  the  sensor  trajectory  policy 
(LO  or  BF),  and  on  the  cycle  length  u.  The  measure  of  effectiveness  (MOE)  used  to  evaluate 
the  LO  and  BF  policies  is  the  long  term  target  detection  rate,  Nd{f)/t  as  t  >  cx). 

3  The  Leap-to-Origin  Trajectory 

In  this  section  we  study  the  leap-to-origin  (LO)  trajectory,  where  the  sensor  instantaneously 
leaps  from  the  destination  to  the  origin  after  completing  a  cycle.  Our  goal  is  to  hnd  the  two 
endpoints  for  the  sensor  to  maximize  the  long-run  detection  rate.  Recall  that  we  assume 
that  Na,  Xi,  and  Ri  are  independent,  that  the  Xfs  have  support  on  3?,  and  that  the  Rfs 
have  support  on  Without  loss  of  generality  (by  changing  the  time  scale  if  necessary) 
we  assume  that  target  arrival  rate  equals  1.  We  also  assume  that  the  sensor  speed  equals  1; 
this  assumption  is  relaxed  later  on  in  Section  5. 

The  problem  of  Ending  the  best  endpoints  is  approached  as  follows: 

(i)  Identify  the  best  origin  point  as  a  function  of  cycle  length. 

(ii)  Find  the  optimal  cycle  length. 

With  regard  to  (i),  define 

H{u)  :=  max{Fx{u-\- a)  -  Fx{a)},  (1) 
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and  let  a*  be  a  maximizer  in  (1).  Observe  that  for  u  arbitrary  and  finite,  Fx{u+a)-Fx{a) 

0  as  |a|  — >  oo,  so  that  — cx)  <  a*  <  oo.  For  a  cycle  length  u,  a*  is  the  optimal  origin  of  the 
sensor,  a*  +  m  is  the  optimal  destination,  and  H{u)  is  the  probability  that  a  target  crosses 
the  perimeter  on  the  sector  {a*,  a*  +  n);  notice  that  H{u)  is  a  distribntion  fnnction  when 
viewed  as  a  fnnction  of  u.  Moreover,  the  hrst  part  of  Assnmption  1  gnarantees  that  H{u)  is 
sufficiently  smooth  for  our  models  to  be  well  posed. 

Lemma  1  Under  Assumption  1,  H{u)  is  continuously  differentiable  and  concave  on  its  sup¬ 
port. 

Lemma  1  follows  almost  immediately  from  the  next  result. 

Lemma  2  Suppose  Assumption  1  holds.  Then  if  fx  has  infinite  support,  the  function  a  : 
3ft  3?  defined  as  the  solution  of 


fx{a{u))  =  fx{a{u)  +  u) 

is  differentiable  and  a{u)  maximizes  Fxiu  +  a)  —  Fx{a). 

Proof  of  Lemma  2.  The  implicit  function  theorem  shows  that  a{u)  is  well  dehned  and 
differentiable.  The  hrst  order  necessary  conditions  for  unconstrained  optimization  show  that 
n  (a  +  u)  -n(a)  =  0  fx{ci)  =  fx{ci  +  u)  and,  since  fx  is  unimodal,  there  is  only 

value  of  a  that  makes  /x(a)  equal  to  fx{ci  +  u)  for  u  given.  0 

Proof  of  Lemma  1.  If  fx  has  inhnite  support,  then 

H'{u)  =  fx{a{u)  +  u){a'{u)  +  1)  -  fx{a{u))a'{u)  =  fx{a{u)  +  u), 

by  Lemma  2.  Therefore,  H'  is  continuous  and  non-increasing,  and  the  lemma  is  proven  for 
the  inhnite  support  situation.  Suppose  now  that  fx  has  hnite  support  from  the  left;  that  is 
i  =  sup{m  ;  Fxiu)  =  0}  is  hnite.  Dehne  as  the  cycle  length  where  the  optimal  trajectory 
origin  hrst  hits  ^.  Therefore,  for  0  <  n  <  we  have  already  shown  that  H{u)  is  continuously 
diherentiable  and  concave;  whereas  for  u  >  ug  we  have  H{u)  =  Fx{i  +  u)  —  Fx{i),  which 
is  easily  seen  to  be  concave  and  continuously  diherentiable.  The  cases  where  fx  has  hnite 
support  from  the  right,  and  from  the  left  and  the  right,  are  similarly  treated.  ® 

A  target  is  detected  on  the  k  +  Vih.  cycle,  for  k  >  1,  when  all  of  the  following  three 
conditions  are  met:  First,  the  target  arrives  to  the  patrolled  sector  (a*  <  Xi  <  a*  +  u)-, 
second,  its  arrival  time  is  at  least  Xi  —  a*  +  {k  —  l)u  (since  we  assume  that  sensor  velocity  is 
1)  and  no  more  than  X^  —  a*  +  ku;  last,  its  reneging  time  is  at  least  Xi  —  a*  ku  —  Ti.  This 
situation  is  illustrated  in  Figure  2,  for  k  =  2:  The  hrst  two  conditions  determine  the  dotted 
region  where  detection  is  possible;  the  last  condition  ensures  that  a  target  that  falls  in  the 
dotted  area  has  a  reneging  time  sufficiently  large  to  be  detected. 

Let  X*  :=  Xi  —  a*  be  the  target  arrival  location  translated  to  a*,  for  i  =  1,2,....  Hence, 
for  an  LO  trajectory  with  origin  a*  and  destination  a*  +  u,  the  expression 

OO 

Nd{{k  +  l)u)  =  Nd{ku)  +  ^I{0  <  X*  <  u-X*  +  {k-l)u  <Ti<  X*  +  ku-Ri  >  X*  +  ku-Ti) 

i=l 
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Figure  2;  Detection  region  for  the  third  cycle  of  an  LO  trajectory,  for  cycle  length  u,  and 
optimal  origin  of  the  sensor  a*. 


represents  the  number  of  detections  at  the  end  of  the  A;  +  I’th  cycle,  for  A;  >  1.  Observe  that 
since  the  hrst  cycle  is  irregular,  we  have  that 

OO 

N,{u)  =  5^/(0  <  X*  <u-T,<  X:-R,  >  X*  -  R) 

i=l 

is  the  number  of  detections  at  the  end  of  the  hrst  cycle. 

We  wish  to  choose  a  cycle  length  and  an  origin  that  maximize  the  detection  rate  Nd{t)/t, 
as  t  >  OO.  The  following  result.  Proposition  4.4.1  in  [14],  provides  theoretical  support  for 
the  computation  of  this  rate. 

Proposition  1  Suppose  that  (Yi  :  i  >  1)  are  random  elements  of  an  Euclidean  space  Ei 
such  that  iV(-)  =  ^  ■)  is  a  Poisson  process  with  Poisson  measure  =  EN{-). 

Suppose  {Ji  :  i  >  1)  are  random  elements  of  a  second  Euclidean  space  E2  with  common 
probability  distribution  F,  and  suppose  N  and  the  p ’s  are  defined  on  the  same  probability 
space  and  are  independent.  Then  the  point  process  Ji)  ^  ')  is  Poisson  with  Poisson 

measure  p,  x  F. 


Proposition  1  guarantees  that  the  process  is  Poisson.  The  property  of  the  Poisson  process 
that  disjoint  intervals  are  independent  implies  that  the  random  variables  that  describe  the 
number  of  detections  in  a  regular  cycle  are  independent  and  identically  distributed,  with 
mean  given  by 


ENd{{k  +  l)u)  -  ENdiku) 

I{0  <  X*  <  u;x*  +  {k  —  l)u  <t<x*  +  ku-,r  >  x*  +  ku  —  t)dFR{r)  dtdFx*{x*) 
Fr{x*  +  ku  —  t)dt  dFx*  (x*) 


'5R  J?R+  J?R+ 

pu  px*-\-ku 


'0  Jx*^{k—l)u 

/u 

FR{t)dt, 


for  A;  >  1,  where  Fr{-)  =  1  —  Fr{-),  and  Fx*{x*)  =  P{X*  <  x*).  Observe  that  although  Nd 
does  not  necessarily  have  a  constant  rate  within  a  cycle,  the  rate  is  periodic  for  all  cycles 
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k  >  2,  and  hence  we  end  up  with  the  above  result.  For  the  hrst  irregular  cycle  we  have 
ENd{u)  =  /  /  /  I{0  <  X*  <  u-,0  <  t  <  x*]r  >  X*  —  t)dFR{r)  dt  dFx*{x*)  <  oo.  (2) 

Equation  (2)  and  the  monotone  convergence  Theorem  imply  that  Nd{u)/t  0  almost 
surely  as  t  >  cx).  Also,  the  Law  of  Large  Numbers  for  Renewal  Reward  processes,  Proposition 
3.4.1  in  [14],  implies  that 

Nd{t)  -  Nd{u)  t-u  FNdijk  +  1)^)  -  FNdjku) 
t  —  u  t  u 

almost  surely  as  t  ^  oo.  Putting  this  all  together, 

Ndjt)  _  Ndju)  ^  Ndjt)  -  Ndju)  t  -  u 
t  t  t  —  u  t 

ENd{{k  +  1)m)  -  ENdiku) 
u 

almost  surely  as  f  >  cx),  for  any  k  >  1.  Therefore,  the  detection  rate  is 

^  ^  ur  JoFR{t)dt 

9lo{u)  :=  H{u) - .  (3) 

u 

Accordingly,  our  goal  in  the  LO  setting  is  to  hnd  a  cycle  length  that  maximizes  gio-  Just 
like  in  (1)  when  selecting  the  best  origin  a{u),  there  may  exist  more  than  one  maximizer  of 
qlo  iu  (3),  in  which  case  we  choose  the  one  with  the  smallest  cycle  length.  That  is,  optimal 
cycle  length  is 

u*  :=  min{M  >  0  ;  gLo{u)  >  gLo{w)yw  >  0}. 

The  term  FR(t)dt  in  (3)  is  the  expected  number  of  detections  in  any  inspected  location 
in  a  cycle  of  length  u,  which  is  then  inflated  by  H{u)  to  account  for  all  surveyed  locations; 
the  u  term  in  the  denominator  is  the  expected  number  of  arrivals  along  a  cycle  of  length  u. 

Remark  1  Changing  the  arrival  rate  of  Na  from  1  to  a  >  0  ehanges  gio  by  «  faetor  a,  but 
does  not  ehange  a*  nor  u* . 

We  dehne  gLoif^)  =  0,  so  that  gio  is  continuous  at  0.  Observe  that  Assumption  1  implies 
gLO  is  continuously  differentiable  in  To  complete  the  notation,  let  f  :=  inf{z/  ;  Rr(zz)  = 
1},  and  u  :=  inf{z/  :  H{v)  =  1}  (both  f  and  u  can  possibly  be  inhnite);  and  write  =  ER 
and  nx  =  EX.  Given  a  differentiable  real  function  i{u)  we  write  i'{u)  as  short  for  di{u)/du. 

We  study  properties  of  the  detection  rate  and  the  optimal  cycle  length  in  various  scenarios 
of  the  LO  trajectory.  First  we  study  several  instances  of  reneging  behavior;  R  bounded  with 
probability  one,  R  deterministic,  and  R  uniformly  distributed.  Second,  we  analyze  the  case 
of  X  bounded,  and  X  uniform.  Last,  we  deal  with  the  general  case  where  both  X  and  R 
have  densities  with  inhnite  support.  We  start  by  analyzing  the  case  where  the  reneging  times 
are  bounded  with  probability  1. 
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Example  1  Suppose  0  <  r  <  oo.  Then 


/  N  _  f  H{u)^°  for  0  <  u  <  r 

'  H{u)^,  foru>f. 


(4) 


Evidently,  if  qlo  is  non-increasing  on  [f,  cx)),  we  must  have  u*  <  f. 

Lemma  3  Suppose  that  0  <  f  <  oo.  Then,  u*  <  f. 

Proof:  Lemma  1  implies  uH'{u)  <  H{u)  for  u  >  f,  so  that  {H{u)/uy  <  0,  and  gfoi'^)  —  0- 


Clearly,  H  strictly  concave  implies  uH'{u)  <  H{u).  The  borderline  case  is  uH'{u)  =  H{u)  for 
u  >  r,  that  is,  H'  constant  on  m  >  f ;  see  Example  4  below.  What  happens  if  uH'{u)  >  H{u) 
for  some  u  >  fl  Then  it  could  happen  that  u*  >  f.  In  any  case,  because  H{u)/u  — 0  as 
M  >  CX),  we  know  that  u*  <  oo. 

Under  the  conditions  of  Lemma  3  it  is  never  benehcial  to  scan  a  sector  longer  than  f 
because  the  detections  gained  by  covering  more  area  would  not  outweight  the  extra  missed 
detections  due  to  reneging.  Two  special  cases  of  Example  1  are  presented  next. 


Example  2  Deterministic  reneging  times:  Suppose  P{R  = 
Equation  (4)  becomes 


groin) 


Hiu),  for  0  <u  <r 
for  u>  f. 


r) 


1,  0  <  r  <  oo. 


Since  gio  non- decreasing  onO  <u<r,  in  view  of  Lemma  3  we  have  u*  =  r. 


Then 


Example  3  If  Ri  f/(0,r),  Equation  (4)  becomes 


arn(u)  =  I  for0<u<f 

[  R(u)£,  foru>f. 


(5) 


Therefore,  u*  <  r,  where  u*  is  the  smallest  maximizer  of  the  top  branch  in  the  RHS  of  (5). 


Our  next  lemma  considers  the  case  of  u  hnite. 


Lemma  4  Suppose  that  0  <  u  <  oo.  Then  u*  <u. 

Proof:  The  fact  that  Fiiit)dt/u  is  non-increasing  in  u,  and  that  Hiu)  =  1  for  all  u  >  u, 
implies  that  u*  <  u.  ® 

Lemma  4  describes  the  scenario  corresponding  to  Xj  having  its  mass  concentrated  in 
some  hnite  interval.  It  says  that  it  is  not  benehcial  to  patrol  further  than  u,  where  no 
targets  arrive.  Our  next  example  illustrates  Lemma  4,  and  shows  that  the  hypothesis  of 
Lemma  3  cannot  be  weakened. 


Example  4  Suppose  Xi  ~  17(0,  m).  Then  Equation  (3)  becomes 


groin) 


HFR(t)dt 

u 


for  b  <u  <u 
,  for  u  >  u. 


Clearly,  u*  =  u.  In  particular,  if  r  <  u,  then  any  cycle  length  in  [r,u]  maximizes  gio- 
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We  now  discuss  properties  of  giaiu)  and  u*  when  both  X  and  R  have  densities  with 
inhnite  support.  First,  it  is  immediate  from  the  dehnition  of  gio  that  gLo{u)u  hr  as 
M  — cx),  so  that  gLoiu)  decreases  towards  0  at  rate  1/m;  see  the  right  part  of  Figure  5. 
Concerning  the  behavior  of  gio  locally  around  0  we  observe  that  gio  is  sub-additive  on 
some  interval  (0,7],  7  >  0: 

fi'Lo(a  +  P)  <  fi'Lo(a)  +  5'lo(/5),  (6) 

for  a,  >  0  and  a  -|-  /?  <  7.  Sub-additivity  follows  because 
(i)  H{a  +  13)  <  H{a)  +  H{{3).  This  follows  since 

H{a  +  (3)  =  max{Fx(a  +  (3  +  a)  —  Fx{a)} 

a 

=  ma,x{Fx{a  +  P  +  a)  —  Fx{a  -|-  a)  -|-  Fx{a  +  a)  —  Fx{a)} 

a 

<  ma,x{Fx{a  +  P  +  a)  —  Fx{a  +  a)}  -|-  max{Fx(Q;  -l-  a)  —  Fx{a)} 

a  a 

=  max{Fx{P  +  a)  —  Fx{a)}  +  max{Fx(a  +  a)  —  Fx{a)} 

a  a 

=  H{p)  +  H{a). 


(ii)  FR{t)dt/{a+P)  <  FR{t)dt/a,  FR{t)dt/ P}.  This  follows  because  FRppdt/u 

is  non-increasing  in  u,  as  can  be  immediately  seen  by  evaluating  its  hrst  derivative. 

Applying  properties  (i)  and  (ii)  leads  immediately  to  (6): 


gLo\oi  +  P)  = - —5 - 

Q  H“  p 


<  {H{a)  -|-  H{P))  min 


^;FR{t)dt  J,^Fn(t)dt 


^  fpFxmt  ,  FRit)dt 
ol  P 

=  fi'Lo(«)  +  gLoiP)- 


In  addition,  gio  is  continuous  at  0.  Therefore,  Lemma  3.6.4  of  [4]  applies,  and  we  have 


gioju)  ^  ^ 
u  ’ 

as  M  — 0,  for  some  A  >  0.  That  is,  g^o  is  locally  linear  around  0;  see  the  left  part  of  Figure 
5.  Thus,  the  behavior  of  gioiF)  for  u  around  0  and  for  all  u  sufficiently  large  implies  that 
the  optimal  cycle  length  must  be  a  stationary  point,  i.e.:  u*  is  the  root  of  the  equation 
g'loi'^)  —  0  which  gio  is  largest.  The  root  problem  can  be  re-written  as 

p'U  pu 

uH'{u)  /  FR{t)dt  =  H{u)  /  FR{t)dt-uH{u)FR{u).  (7) 

Jo  Jo 
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The  last  equation  says  that  we  increase  cycle  length  until  the  marginal  increase  in  target 
detections  equals  the  marginal  decrease  in  lost  target  detections  due  to  reneging.  In  general, 
however,  it  is  not  possible  to  analytically  solve  Equation  (7),  so  standard  numerical  methods 
must  be  utilized;  the  computational  cost  grows  with  the  number  of  roots  that  need  to  be 
found  in  (7). 

Example  5  Suppose  Xi  ~  iV(0, 1)  and  Ri  ~  exp(l).  Then  a*  =  —1.025,  u*  =  2.05  is  the 
root  of  g'l^Qiu)  =  0,  and  groir*)  =  -295;  see  Figure  5. 

Remark  2  Concerning  Example  5,  an  exponential  reneging  distribution  leads  to  the  same 
solution  that  arises  in  the  following  setting:  Upon  arrival,  each  target  emits  a  signal  of 
constant  amplitude  that  decays  exponentially  in  time  with  rate  9;  this  is  known  as  a  shot 
noise  process  [If]-  This  would  be  the  case,  for  example,  if  the  footprint  and  marks  left  by 
the  infiltrators  disappear  at  an  exponential  rate  in  time.  In  particular,  the  signal  strength  of 
arrivals  detected  in  the  {k  +  1)  ’th  cycle  is  given  by 

OO 

<Xi<a*  +u,Xi  +  {k-  l)u  <  +  ku). 

i=l 

Following  the  approach  taken  in  [If],  P-  326,  we  find  the  expected  average  signal  strength 
of  targets  detected  in  an  arbitrary  cycle  is  H{u)  J]]  exp{—9t)dt/u,  which  is  gLoiu)  when 
Ri  ~  exp(6*).  In  other  words,  a  system  where  targets  that  arrive  according  to  a  shot  noise 
process,  has  the  same  solution  as  the  system  where  targets  renege  according  to  an  exponential 
distribution. 

Next,  we  study  the  effect  of  variability  in  the  targets’  arrival  location.  The  question  we 
have  in  mind  is:  Does  an  increase  in  :=  var(X)  cause  g^o  fo  decrease?  The  answer  is 
not  necessarily,  as  can  be  seen  in  the  following  example. 

Example  6  Suppose  that  X  has  a  density  that  is  a  mixture  of  a  uniform  and  a  Pareto 
density, 

f  {  \  _  j  1)  for  0  <  X  <9 

forx>9, 

and  that  P{R  <  9)  =  1,  for  0  <  9  <  1.  It  is  easy  to  see  that  II{u)  =  l  +  —  9)/u  for  u  >  9, 

and  that  9  >  1/2  implies  {HipTj/u)'  <  0  for  u  >  9.  Therefore,  according  to  Lemma  3  and 
Example  f  we  have  u*  =  9,  and  gLoiu*)  =  Tr-  Compare  this  with  the  case  X*  ~  U(0,u), 
where  u  >  9.  In  this  case  we  have,  according  to  Example  f,u*  =  9  and  gLoiu*)  =  Tr/H. 
Therefore,  ifu>  1,  the  optimal  detection  rate  for  the  mixture  density  case  is  larger  than  that 
of  the  non-uniform  case.  However,  a\  in  the  mixture  case  is  infinite  (due  the  heavy-tails 
of  the  Pareto  density),  while  in  the  second  case  is  finite.  This  example  lets  us  conclude 
that  larger  arrival  location  variability  does  not  necessarily  lead  to  lower  detection  rate  for  the 
optimal  cycle  length. 

However,  for  re-scaled  versions  of  X  it  is  true  that  lower  variance  leads  to  greater  detection 
rate.  To  show  this  it  suffices  to  compare  the  respective  H{-)  terms. 
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Lemma  5  Let  X  =  Xjox,  a,*  a  maximizer  of  P{a  <  X  <  a  +  u),  and  H{u)  =  P{d*  <  X  < 
d*  +u).  Suppose  0  <  ax  <  oo.  Then  H{u)  <  H{u)  if  and  only  if  ax  hi¬ 
proof:  For  (Tx  >  1  we  have  H{u)  =  P{axd*  <  X  <  axd*  +  axu)  >  P{a*  <  X  < 

a*  +  axu)  >  P{a*  <  X  <  a*  +  u)  =  H{u).  On  the  other  hand,  for  ax  <  1  we  have 
H{u)  <  P{axd*  <  X  <  axd*  +  u)  <  P{a*  <  X  <  a*  +  u)  =  H{u).  ® 

Example  7  Suppose  Xi  ~  iV(0,a^),  and  Ri  ~  exp(l).  The  top  of  Figure  3  illustrates  the 
optimal  detection  rate  as  a  function  of  a;  as  expected  from  Lemma  5,  gLoip*)  decreases  as 
a  increases.  The  bottom  of  Figure  3  shows  the  optimal  cycle  length  increasing  with  a.  Not 
shown  is  the  optimal  origin  a*  =  —u*/2. 


Figure  3:  Detection  rate  and  cycle  length  as  a  function  of  a. 


4  The  Back- and- Forth  Trajectory 

In  this  section  we  study  the  back-and-forth  (BF)  trajectory — a  scenario  that  arises  when 
sensors  are  mounted  on  platforms,  such  as  ground  patrol  units  or  UAVs.  Our  goal  is  to 
analyze  cycle  trajectories  and  hnd  an  optimal  cycle  length  in  order  to  maximize  the  long-run 
detection  rate  for  targets  before  they  renege. 

The  setting  is  the  same  as  that  in  the  LO  scenario:  X^,  R^,  and  Na{-)  are  mutually 
independent;  the  Xj’s  and  Rfs  are  continuous  random  variables  whose  distributions  satisfy 
Assumption  1;  and  the  Rfs  take  values  on  3?+.  The  speed  of  the  sensor  is  1,  and  the  targets 
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Figure  4:  Detection  region  for  the  first  regular  cycle  of  a  BF  trajectory,  for  cycle  length  2u, 
and  origin  of  the  sensor  a. 


arrive  according  to  a  Poisson  process  with  rate  1.  As  shown  in  Figure  1,  because  the  re-visit 
time  to  the  center  of  the  trajectory  is  shorter  than  that  of  the  endpoints,  targets  renege  more 
often  in  an  area  closer  to  the  trajectory  endpoints. 

In  order  to  obtain  an  expression  for  the  long-run  detection  rate,  we  take  a  similar  approach 
to  the  one  taken  to  determine  qlo-  Consider  Figure  4:  We  hlled  with  ”  marks  the  region 
where  targets  get  detected  (if  no  reneging  occurs)  when  the  sensor  moves  from  a  +  u  to  u, 
and  we  hlled  with  marks  the  region  where  targets  get  detected  (if  no  reneging  occurs) 
when  the  sensor  moves  from  u  back  to  a  -|-  u.  The  cycle  length  is  2u,  and  a  is  the  cycle 
origin.  For  a  target  i  to  fall  in  the  ”  marked  region  and  get  detected  it  must  satisfy  three 
conditions:  a  <  Xi  <  a  +  u,  Xi  —  a  <  Ti  <  2u  —  (X*  —  a),  and  Ri  >  2u  —  Ti  —  [Xi  —  a). 
On  the  other  hand,  for  a  target  to  fall  in  the  marked  region  and  get  detected,  it  must 
satisfy:  a  <  X^  <  a  +  u,  2u  —  {Xi  —  a)  <  R  <  2u  +  (Xj  —  a),  and  Ri  >  2u  —  R  +  {Xi  —  a). 
A  similar  analysis  can  be  carried  out  for  all  the  regular  cycles,  by  introducing  a  cycle-index 
k,  just  like  we  did  in  the  LO  case.  Hence,  Proposition  1  ensures  that  the  detection  process 
for  the  BF  trajectory  is  Poisson,  so  that  the  number  of  detections  in  each  regular  cycle  is 
a  sequence  of  independent  and  identically  distributed  random  variables,  with  the  mean  for 
each  cycle  given  by 


f*2u—{x—a) 


Fr{2u- 


(X 


a)—t)dt  dFx{x)  -|- 


pa-\-u  p2u+{x—a) 


'  2u—(x—a) 


Fji{2u  —  t+{x  —  a))dt  dFx{x) 


pa-\-u  p2{u—(x—a)) 


FR{t)dt  dFx{x)  + 


r*a+4i  p2[x—a) 


FR{t)dt  dFx{x) 


'a  v/  0 

/u 

fx{x  +  a)c{x,  u)dx, 


'a  v/  0 


where  fx  is  the  density  of  the  Xj’s  and 


p2{u—x)  p2x 

c{x,u)  ■.=  /  FR{t)dt+  /  FR{t)dt,  (8) 

Jo  Jo 

ior  0  <  X  <  u. 

Since  the  mean  number  of  detections  in  the  hrst  irregular  cycle  is  hnite,  the  Law  of  Large 
Numbers  for  Renewal  Reward  processes  ensures  that  the  long-run  detection  rate  for  targets 
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is  given  by 


9bf{u) 


maxa{  fx{.x  +  a)c(x,  u)dx} 
2u 


(9) 


The  constant  (i*bf  ~  a{u)  maximizes  fxix  +  a)c{x,  u)dx,  and  is  the  best  trajectory  origin 
for  cycle  length  u.  Since  f(x  +  a)  — 0  as  |a|  — >  cx)  and  c{x,  u)  <  2/i/j  for  0  <  a;  <  we  must 
have  |a*|  <  cx).  From  Equation  (8),  it  can  be  seen  that  the  function  c{-,u)  is  symmetric  for 
each  M  >  0  about  x  =  u/2.  In  addition,  differentiation  with  respect  to  x  shows  that  c(-,  u)  is 
increasing  and  concave  in  x  for  hxed  u.  Dehne  gsFib)  =  0  so  that  qbf  is  continuous  in  3?+, 
and  observe  that  Assumption  1  implies  that  Qbf  continuously  differentiable  in  3?+  (the 
proof  of  this  fact  mimics  that  of  Lemma  1).  The  optimal  cycle  length  is  the  smallest  cycle 
length  that  maximizes  Qbf'- 


u*BF  ■=  niin{M  >  0  ;  QBFiu)  >  gBF{w),'iw  >  0}, 

The  analysis  in  this  section  is  similar  to  that  of  the  LO  trajectory,  and  we  draw  compar¬ 
isons  between  LO  and  BF  whenever  appropriate.  In  this  regard,  the  next  lemma  shows  that 
the  BF  trajectory  can  do  no  better  than  the  LO  trajectory,  for  any  cycle  length  u. 


Lemma  6  For  all  u  >0,  qbf^u)  <  gLoiu). 


Proof:  Observe  that  Fr  non-increasing  implies  that  c{x,u)  <  2  FB(t)dt,  for  all  0  <  a;  < 
u.  Hence 


dEF^u)  — 


maXa{C  fx{x  +  a)c{x,  u)dx}  maXa{C  fx{x  -f  a)dx}  C  FR(t)dt 


2u 


< 


u 


=  9lo{u). 


In  addition. 

Corollary  1  Lemma  6  implies 


9bf{ubf)  ^  9lo{ubf)  <  9lo{ufq), 


that  is,  the  optimal  detection  rate  for  the  LO  trajectory  cannot  be  smaller  than  the  optimal 
detection  rate  for  the  BF  trajectory. 


We  now  focus  on  properties  of  the  optimal  cycle  length  u*bf-  We  start  with  the  analogue 
of  Lemma  3,  dealing  with  the  case  of  Ri  having  hnite  support,  P{Ri  <  f)  =  1,  for  0  <  f  <  cx. 
The  key  idea  is  that  for  u  3>  f,  gBriu)  ^  9 lo{u).  Specifically,  for  m  >  f  we  have 


„  ,  Fx{u  +  a*BF)  -  Fx{a*BF) -f{u) 

9bf[u)  —  Hr - 

u 

(cf.  Equation  (4))  where 


(10) 


ip{u)  :  = 

fx{x  +  a*RF)  (^Hr  -  lo""  FR{t)dt^  dx  +  fx{x  +  (^hr  -  FR{t)dt^  dx 

2/i_R 


is  the  “error”  term. 
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Lemma  7  Suppose  that  0  <  r  <  oo,  that  qbf  is  differentiable  in  and  that  the  error 
term  (p{u)  satisfies 

u  ifxiu  +  Obf)  -  <  Fxiu  +  a*BF)  -  -  ip{u),  (11) 

for  all  u  >  f.  Then  u*Bp  <  f. 

Proof:  The  hypothesis  implies  that  —  0  m  >  f.  ® 

Remark  3  Condition  (11)  is  the  same  as  that  of  Lemma  3,  but  for  the  error  term  ip.  The 
error  term  becomes  insignificant  as  f  becomes  large,  but  may  be  relevant  for  small  values  of 
f — see  Example  (8)  below. 

Example  8  When  P{R  =  f)  =  1  for  some  constant  f  >  0  we  have 

c{x,  u)  =  min{2a;,  f}  +  min{2(M  —  x),f}. 

For  u  <  fj^  we  have  c{x,u)  =  2u,  and  gBFiu)  =  Fxiu  +  afip)  —  Exiafip)  non- decreasing 
implies  f/2  <  u*bf-  U  the  conditions  of  Lemma  7  apply,  we  have  f/2  <  u*bf  <  L . 

The  analogue  of  Lemma  4  to  the  BF  setting  is  straightforward. 

Lemma  8  Suppose  P(n  <  X  <  it)  =  1,  where  u  =  sup{M  :  P(X  >  m)  =  1}  and  it  =  inf{M  : 
P(X  <u)  =  1},  and  —oo  <  u<  it  <  oo.  Then  u*bf  <u  —  u. 

Proof:  For  u^>u  and  <  u,  we  have 

lu^  fxix)cix,u)dx  ^  fxix)cix,u)dx  ^  £  fxix)cix,u)dx 
2(m+  — M_)  2(m+  — M_)  “  2(-u  —  m) 

The  other  possibilities, 


•  Up  —  U-  >  u  —  u,  U-  >  u,  and  >  u, 


•  u+  —  U-  >  u  —  u,  U-  <  u,  and  u+  <  u, 


are  treated  similarly.  ® 

Regarding  the  general  case  when  both  Xi  and  Ri  have  inhnite  support,  the  results  in  the 
BF  case  match  those  of  the  LO  case: 

(i)  From  the  dehnition  of  qbf  we  have  qbf  decreasing  at  rate  as  m  — >  cx).  This 

implies  that  graiu)  —  gBriu)  — 0  as  m  — >  cx);  see  the  right  part  of  Figure  5. 

(ii)  gBF  is  sub-additive;  the  proof  of  this  fact  is  similar  to  the  LO  one.  Hence,  gBF  is  locally 
linear  (increasing)  around  0.  Moreover,  it  can  be  shown  that  the  linear  constants  are 
the  same  in  both  the  LO  and  BF  cases. 

(iii)  Items  (i)  and  (ii)  imply  that  the  optimal  cycle  length  satishes  g^Fi"^)  —  0)  which  can 
be  solved  numerically  at  a  cost  that  grows  with  the  number  of  roots. 
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(iv)  Effects  of  arrival  location  variability:  Here  we  also  have  that  re-scalings  of  X,  such  as 
cX,  lead  to  lower  detection  rate  if  and  only  if  the  scaling  constant  c  is  greater  than  1. 

The  next  example  illustrates  many  of  the  results  of  this  section. 

Example  9  Suppose  Xj  ~  X(0, 1)  and  Ri  ~  exp(p),  for  p  >  0.  We  have  =  —u/2  and 
completing  squares  in  equations  (9)  and  (8)  results  in 

,  ,  ($(m/2)  -  $(-m/2))  -  e^P^-P^(^(2p  +  u/2)  -  $(2p  -  u/2)) 

gBF[u)  = - . 

up 

Figure  5  shows  the  detection  rates  for  the  LO  case  on  top  and  for  the  BF  case  on  bottom 
when  p  =  1.  As  expected  from  Lemma  6,  we  observe  gspiu)  <  gLoiu)  for  all  u  >  0.  In 
addition,  gspiu)  ~  gLoiu)  locally  around  0  because  there  are  no  reneges  there,  and  gspiu)  ~ 
QLoidj)  for  all  u  sufficiently  large  because  the  expected  number  of  detections  per  cycle  converges 
to  pr  (i.e.:  the  numerator  in  Eqs.  (3)  and  (9)  converges  to  pn  as  u  ^  oo),  yet  the  increasing 
cycle  length  causes  the  detection  rate  to  decrease  at  speed  0{u~^).  The  operational  take-away 
is  that  there  is  not  a  significant  difference  between  the  LO  and  BF  trajectories  when  the  cycle 
length  is  small  or  large. 


Figure  5:  Detection  rate  for  LO  and  BF  trajectories. 

Figure  6  presents  u^p,  gBp{u*Rp),  u*pQ,  and  gpoiFfo)  as  a  function  of  the  reneging  rate 
p.  We  observe  in  the  top  part  that  Upp  >  u*pQ  for  all  p;  the  reason  for  this  is  that  the  BF 
trajectory  compensates  the  reneges  at  the  endpoints  with  the  extra  detections  provided  by  a 
larger  trajectory  length.  (This  argument  can  be  made  rigorous,  but  there  is  not  much  value  in 
comparing  optimal  cycle  lengths).  Regarding  the  bottom  part  of  Figure  6  we  have,  as  expected 
from  Corollary  1,  that  gBp{u*pp)  <  gLoixfo),  with  the  difference  gLoi^fp)  —  gBp{u*pp)  being 
largest  for  moderate  values  of  p.  In  other  words,  when  targets  have  an  expected  renege  time 
that  is  sufficiently  small  or  large,  the  LO  trajectory  is  not  much  better  than  the  BF  trajectory. 
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Figure  6:  Cycle  length  and  detection  rate  as  a  function  of  reneging  rate. 

5  Effects  of  Velocity  and  Investigation  Time 

Although  the  baseline  scenarios  described  in  the  previous  sections  are  idealized,  they  can  be 
adjusted  in  ways  that  better  approximate  sensor  operation  in  the  held.  In  this  section  we 
put  forward  two  such  extensions.  In  the  hrst  extension  we  study  the  effect  of  sensor  velocity 
on  the  long  term  detection  rate.  In  particular,  we  consider  the  tradeoff  between  visiting  a 
point  on  the  perimeter  more  frequently,  which  mitigates  the  effect  of  reneging,  and  missing 
targets  due  to  decreased  detection  capability  when  sensor’s  velocity  increases.  In  the  second 
extension  we  study  the  case  where  it  takes  time  to  investigate  each  detected  target. 

Remark  4  The  effects  of  these  extensions  for  both  the  LO  and  BF  settings  are  qualitatively 
similar.  In  this  paper  we  only  present  the  LO  trajectory  since  it  is  more  tractable. 

5.1  Effects  of  velocity 

If  sensor  velocity  has  no  impact  on  its  ability  to  detect  inhltrators  that  cross  its  trajectory 
then  it  is  obvious  that  the  optimal  cycle  length  and  detection  rate  increase  with  sensor 
velocity,  and  maximum  velocity  is  optimal.  On  the  other  hand,  if  high  velocity  leads  to  a 
lower  detection  capability,  low  sensor  velocity  may  be  optimal.  In  this  section  we  investigate 
the  effect  of  this  dependency. 

Given  that  the  sensor  looks  at  a  target,  we  consider  the  case  where  the  target  is  detected 
with  a  probability  that  is  non-increasing  in  the  sensor  velocity;  this  case  captures  the  situa¬ 
tion  where  the  sensor  operator  tends  to  miss  targets  as  the  sensor  velocity  increases.  Because 
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the  detection  process  in  this  setting  is  a  thinned  version  of  the  detection  process  iVrf(-)  (see 
[14],  p.  317),  in  order  to  maximize  detection  rate,  velocity  v  and  cycle  length  u  should  be 
selected  in  the  LO  context  so  as  to  maximize 


9lo{u,v) 


r]{v)H{u)  FR{t)dt 
u/v  ’ 


(12) 


where  r]{v)  is  the  detection  probability  when  the  sensor  moves  at  velocity  v.  The  thinning  of 
the  detection  process  shows  as  the  r]{v)  term  in  Equation  (12).  Naturally,  we  seek  an  answer 
to  the  following  question:  What  is  the  (velocity,  cycle  length)  pair  that  maximizes  detection 
rate?  We  assume  that  the  detection  probability  r]{v)  is  differentiable  for  all  n  >  0. 


Example  10  Suppose  rjiv)  =  p  <1  is  positive  and  eonstant  for  all  v  >  0.  This  is  the  ease 
where  sensor  veloeity  has  no  impact  on  the  sensor  ability  to  detect  targets.  We  have 


9lo{u,v) 


pH{u)  FR(t)dt 
u/v 


For  u  arbitrary,  dgLo{u,v)/dv  >  0  shows  that  going  at  maximum  speed  (say  v)  is  optimal. 
Changing  variables,  ip  =  u/v,  shows  that  maximizing  gLo{u,v)  is  the  same  as  maximizing 

pH{vi))  FR{t)dt 

Hence,  we  are  in  the  framework  of  Section  3,  with  tjj  in  place  of  u.  Let  ifj*  be  the  smallest 
maximizer  in  (13),  and  u*{v)  =  vfj* .  A  repeat  of  the  arguments  used  in  Section  3  shows 
that  u*{v)  >  u*  and  gLoiT*{/v),v)  >  gLoiu*)  if  and  only  if  v  >  1.  In  other  words,  increasing 
velocity  can  only  improve  detection  rate. 


The  analogue  of  Example  1  is 


Example  11  Suppose  0  <  r  <  oo.  Then 


9lo{u,v) 


p{v)H{u)^ — 0  — 
p{v)H{u)^,  for  u/v  >  f. 


In  the  above  example,  if  H  is  strictly  concave  then  H{u)/u  is  decreasing  for  all  m  >  0,  and 
so  we  must  have  u*  <  v*f  (cf.  Lemma  3).  The  continuation  of  Example  2  is 

Example  12  Suppose  that  P{R  =  f)  =  1,  0  <  f  <  oo.  Then,  if  H  is  concave,  we  have 
u*{v)  =  fv  and 

9lo{u*{v),v)  =  p{v)H{f). 

Because  p{v)H{f)  b  if  p{v)  b  as  v  ^  oo,  we  must  have  b  <  v*  <  oo,  so  that  the 
optimal  velocity  satisfies  {p{v)H{vf)y  =  0. 


The  following  is  an  example  of  practical  interest  because  it  treats  the  case  of  a  uniform 
distribution  of  the  targets’  arrival  location  along  the  perimeter. 
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Example  13  Suppose  Xj  ~  U{0,u).  For  any  v  >  0,  Example  4  implies  that  u*{v)  =  u,  so 
that 


V{v)v  Fji{t)dt 


9lo{u,v)  = 


u 


Observe  that  for  all  0  <  v  <  e,  e  >  0,  we  have 

gLoiu,v)  <  e 


u 


which  implies  that  gLo{u,v)  — 0  as  v  — >  0.  Similarly,  for  all  v  >  f  >  0, 


9lo{u,v)  <  - z - =  viO: 

u 

shows  that  gLoiu,v)  — 0  as  v  — >  cx) .  Therefore,  the  optimal  velocity  is  0  <  v*  <  oo,  and  can 
be  found  by  solving  dgLoih,v) /dv  =  0.  In  this  regard,  gLoih,v)  strictly  concave  in  v  would 
imply  the  existence  of  only  one  such  root. 


Looking  now  at  the  general  case,  we  start  by  exploring  some  properties  of  the  optimal 
(velocity,  cycle  length)  pair. 

Lemma  9  Suppose  that  rj{v)  is  non-increasing  and  goes  to  0  as  v  ^  oo.  Then  0  <  u*  <  oo 
and  0  <  V*  <  oo. 


Proof: 


Observe  that  for  all  m,  n  >  0  we  have 


FR(t)dt 

u/v 


<  1. 


(14) 


We  study  four  complementary  scenarios:  For  0  <  m  <  e,  e  >  0,  Equation  (14)  implies 
gLoiT,v)  <  H{e)rj{fS).  Second,  for  0  <  n  <  e,  e  >  0,  we  have  gLo{u,v)  <  eri{0)p,fi/u. 
Third,  for  m  >  .^  >  0:  gLoiT,v)  <  vri{0)p,ji/f.  Last,  for  n  >  .^  >  0,  Equation  (14)  implies 
gLoiu,v)  <  rj{ff).  Sending  e  — 0  and  cx)  in  the  above  scenarios  completes  the  proof.  ® 
The  last  lemma  lets  us  conclude  that  the  optimal  {u* ,v*)  pair  must  be  a  stationary  point 


dgLoiu,v)  dgLoiu,v) 

- - - =  0,  and  - - - 

ou  ov 


(15) 


Just  like  in  the  LO  and  BE  case,  we  use  numerical  optimization  to  hud  such  stationary 
point,  and  select  the  one  for  which  the  detection  rate  is  largest.  Finding  the  optimal  [u*,  v*) 
pair  can  become  numerically  costly  when  there  exist  many  roots  that  need  to  be  found  and 
compared. 


Example  14  Suppose  that  Xi  ~  exp(l),  Ri  ~  1/(0, 1),  and  rjiy)  =  exp(— v).  Then  a  direct 
calculation  shows  that 

g{u,v)  =  e“’^(l  —  e“'“)(l  —  .bu/v), 

for  u  <  v;  this  constraint  is  justified  by  the  discussion  following  Example  11.  Solving  (15) 
yields  {u*,v*)  =  (.610,  .726),  and  g{u* ,v*)  =  .128.  We  now  explain  how  this  .128  fraction  is 
obtained.  P{X  <  u*)  =  .457  is  the  fraction  of  targets  that  arrive  to  the  part  of  the  perimeter 
that  is  searched.  Of  these,  .457  x  FR{t)dt  =  .223  are  observed  before  reneging.  Finally, 
.223  X  exp{—v*)/{u*/v*)  =  .128  is  the  fraction  of  targets  that  is  detected. 
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Figure  7:  The  LO  trajectory  (left)  and  BF  trajectory  (right),  subject  to  investigation  time. 

5.2  Investigation  time 

So  far  we  have  postulated  operational  scenarios  under  the  assumption  that  the  detection 
occurs  instantaneously.  In  certain  situations,  however,  targets  are  investigated  by  the  sensor 
upon  detection,  during  which  time  the  sensor  remains  immobile.  This  type  of  sensor  policy  is 
common  when  the  identihcation  of  the  target  is  reqnired  before  a  decision  is  made  as  to  the 
necessary  reaction,  or  when  the  sensor  has  to  escort  the  target  nntil  a  response  force  arrives 
at  the  scene.  Fignre  7  illustrates  this  scenario  when  the  investigation  time  is  constant. 

As  in  the  previous  sections,  we  would  like  to  hnd  a  cycle  length  that  maximizes  detection 
rate  as  t  ^  oo.  Unfortnnately,  for  a  given  cycle  length  u,  the  detection  rate  cannot  be 
analytically  computed,  and  must  be  estimated  by  Monte  Carlo  simulation.  This  happens 
because  disjoint  intervals  do  not  have  independent  number  of  detections,  so  the  detection 
process  is  not  Poisson  anymore.  Hence,  we  must  use  simulation  optimization  ideas  to  hnd  a 
good  cycle  length. 

To  make  our  developments  simple  we  assnme  that  investigation  times  are  constant  and 
eqnal  to  6^  >  0.  The  key  to  approximating  detection  rate  for  large  valnes  of  t  is  to  notice  that 
the  nnmber  of  detections  that  fall  within  complete  nninterrnpted  trajectories  and  the  times 
elapsed  between  these  trajectories  are  i.i.d.  random  variables.  In  other  words,  the  system 
regenerates  whenever  a  trajectory  with  no  detections  occurs. 

To  get  started,  let  Z  =  {Zi  :  i  >  0)  be  a  stochastic  process  with  state  space  on  . 
A  realization  Zi  =  (ti,o;  •  •  •  >  ^i,n)  says  that  the  Fth  trajectory  starts  at  time  ti^,  and 
the  j’th  detection  occurs  at  time  tjj,  for  j  =  1, . . .  ,n,  n  being  the  last  detection  in  that 
trajectory.  We  set  Zq  =  (0)  to  indicate  that  there  are  no  detections  in  the  hrst  (irregular) 
cycle.  Because  we  know  velocity  and  the  endpoints  of  the  patrolled  sector,  given  Zi  we 
also  know  the  location  of  the  arrivals  in  that  cycle.  Observe  that  the  ZiS  are,  in  general, 
not  i.i.d.  random  variables.  However,  homogeneous  Poisson  arrivals  imply  that  P{Zn  G 
■\Zn-i, . . . ,  Zi)  =  P{Zn  G  ■\Zn-i),  and  hence  the  process  Z  is  Markov.  Let  us  dehne  stopping 
times  Ti  =  inf{i  >  0  :  \Zi\  =  1},  and  =  inf{i  >  r„_i  :  \Zi\  =  1}  for  n  >  2,  where  by  \Zi\ 
we  denote  the  cardinality  of  Zi.  Observe  that  the  r^s  are  i.i.d.  random  variables.  Observe 
that  if  targets  never  renege  it  is  possible  to  have  r  =  oo.  However,  when  0  <  /tr  <  cx),  we 
have  the  following 

Lemma  10  If  0  <  p,R  <  oo  then  Eti  <  oo. 

Proof:  Assume  momentarily  that 

P(ti  >  k)<il-  exp(-fiR))'‘, 


(16) 
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CX). 


for  k  >  1.  Hence, 

OO 

Eti  =  >  k)  <  exp(/iR)  < 

A:=0 

To  complete  the  proof  we  use  induction  to  show  (16): 

P(ri  >  1)  =  P(|Zi|  >  1)  =  1  -  exp  (-H{u)  [  FR{t)dt  1  <  1  -  exp(-p/j), 


where  the  second  equality  follows  from  the  fact  that  the  process  is  Poisson.  A  similar 
argument  verihes  the  inequality 


P{\Zk+i\  >  l\\Zk\  >  1)  <  1  -  exp(-pR),  (17) 

for  k  >  1.  Assume  that  P(ri  >  /c)  <  (1  —  exp(— p/j))^  for  k  >  1.  Then, 

Piji  >  /c  +  1)  =  P{\Zk+i\  >  1;  \Zk\  >  1,  •  ■  ■  ,  \Zi\  >  1) 

=  P{\Zk+i\  >  l||^fc|  >  E)P{\Zk\  >  I,--  -  ,lZil  >  1)  (Markov  property) 

<  (1  —  exp(— /iij))(l  —  exp(— pr))^  (by  (17)  and  the  induction  hypothesis). 


We  will  also  make  use  of  the  following 


Corollary  2  If  0  <  <  oo  then  E\Zk\  <  oo,  for  k  >1. 


Proof:  The  proof  follows  from  (17)  since  P{\Zk\  >  1)  =  E{P{\Zk\  >  l||Zfc_i|  >  1))  < 

1  —  exp(— pr)  implies  E\Zk\  <  oo. 

Lemma  10  and  its  corollary  let  us  use  the  Law  of  Large  Numbers  for  Renewal  Processes 
to  conclude  that 


if  detections  by  time  t 


t 


almost  surely  as  t  >  cx),  where 


giNviu) 


giNviu), 

(18) 

1-1)) 

1))  -f  uEti' 

(19) 

The  term  — 1)  is  the  (random)  number  of  detections  in  a  renewal,  and  6*  — 

1)  +  UTi  is  the  (random)  renewal  time.  Lemma  10  and  its  corollary  ensure  that  giNviu)  is 
well-dehned.  Observe  that  giNv{u)  — 1/9  as  the  target  arrival  rate  goes  to  cx,  for  p,  <  oo. 
Although  the  notation  does  not  reflect  it,  the  right-hand  side  in  (19)  depends  on  cycle  length 
u. 

Assume  that  giNV  has  a  unique  maximizer  uj^y.  Then,  because  giNviu)  <  gioiu)  for 
all  M  >  0,  we  have  giNviu)  =  0{u~^)  and  other  words,  the  search  space 

of  the  optimal  cycle  length  is  bounded.  Regarding  the  optimal  cycle  origin,  when  f  <  u  it 
is  possible  to  show  that  a*  is  a  solution  of  (1).  Unfortunately,  the  problem  is  analytically 
intractable  in  the  general  case  when  f  >  u  (possibly  f  =  cx).  We  carried  out  an  exhaustive 
simulation  optimization  analysis  for  various  distributions  of  R  and  X  that  meet  Assumption 
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1,  and  found  that  a*  is  a  solution  of  (1)  in  all  scenarios.  This  leads  us  to  conjecture  that  the 
optimal  cycle  origin  solves  (1). 

In  order  to  hnd  an  estimate  of  a  cycle  length  that  maximizes  gmv,  we  employ  the  method 
of  Sample  Average  Approximation  (SAA);  see  chapter  6  of  [15]  for  an  overview  of  SAA,  and 
Theorem  2.1  of  [16]  for  the  version  of  SAA  used  in  this  paper.  The  basic  idea  of  SAA  is  to 
maximize  giNv{u)  by  optimizing  (the  ratio  of)  sample  averages 


gniu)  := 


Vn{u) 

WniuY 


where  Vn{u)  =  1/n  X)r=i  ^niu)  =  {Vi{u),Wi{u)ys  are  i.i.d. 

random  samples  drawn  from  the  same  distribution  of  “1); 

Although  this  distribution  is  generally  unknown,  we  can  obtain  samples  from  it  by  generating 
a  sufficiently  large  stream  of  reneging  times  Ri,  locations  Xj,  and  arrival  times  T*,  and  then 
using  discrete  event  simulation  [11]. 

In  order  to  use  SAA,  we  must  have  V (u)  and  W (u)  Lipschitz  continuous: 


|I^(Mi,a;)  -  V{u2,uj)\  <  i'{uj)\ui  -  M2I,  (20) 

and 

\W{ui,uj)  -  W{u2,uj)\  <  ^jJ{uJ)\ul  -  U2I,  (21) 

where  a;  is  a  sample  path,  and  u  and  ip  are  functions  of  uj  with  hnite  second  moments.  It 
is  possible  to  bound  the  LHS  of  (20)  by  the  number  of  arrivals  over  the  largest  regenerative 
cycle  Ti  (amongst  the  two  cycles  stemming  from  selecting  ui  or  U2),  which  due  to  the  Poisson 
arrival  assumption  and  Corollary  2  has  a  hnite  second  moment.  Hence,  for  |ui  —  M2I  >  1, 
Equation  (20)  is  satished.  Using  the  same  type  of  argument,  it  is  possible  to  bound  the 
LHS  of  (21),  for  \ui  —  M2I  >  1-  Unfortunately,  this  approach  does  not  work  for  (20)  and 
(21)  when  \ui  —  U2I  <  1,  and  our  setup  is  too  complex  to  enable  analytical  verihcation  of 
both  Lipschitz  conditions.  To  overcome  this  difficulty,  we  ran  a  simulation  with  a  very  large 
number  of  replications  and  several  choices  of  parameters  Ui  and  U2  such  that  \ui  —  M2 1  <  1; 
the  simulation  results  indicate  both  Lipschitz  conditions  are  satished  when  \ui  —  M2I  <  1- 
Let  iin  be  the  cycle  length  obtained  by  numerically  maximizing  pn-  Then,  Theorem  2.1 
of  [16]  (a  version  of  the  Delta  Theorem)  ensures  that 


n 


1/2  [  9n{Un)  —  giNv{u*jjyiy) 


(Jn 


^iV(0,l), 


as  n  — cx),  where 


=  (Ji, 


I  2  ^  Vn{Un) 


^W^iUn)  ' 

and  =  (n  -  1)“^  E*=i(^*(“n_)  -  K(^^n))^  =  (n  -  1)“^  “  ^n(^^n))^ 

4^^  =  {n-  1)-1  -  %{Un)){WYUn)  -  Wn{Ur,)). 

This  result  provides  theoretical  support  for  building  asymptotically  valid  conhdence  in¬ 
tervals  for  giNviu*iNv)- 

P  {giNvi'^INv)  ^  [Pn:  Rn])  —^1  —  5, 
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Investigation  time  6 


Parameter 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

gn{Un)  ±  1.96(Jn/n^/^ 

.295 

.287  ±  .01 

.276  ±  .008 

.265  ±  .008 

.244  ±  .007 

.240  ±  .007 

Un 

2.05 

2.04 

1.97 

1.90 

1.92 

1.90 

dn 

-1.02 

-1.02 

-.99 

-.95 

-.96 

-.96 

Table  1.  Detection  rate  with  different  investigation  times. 


as  n  — cx),  where 


Ln  —  9n{Un) 


dn 


R-n  fl'n(hn) 


+  ;2 


dn 


and  2:  is  selected  so  that  P{—z  <  N{0, 1)  <  2:)  =  1  —  5. 


Example  15  Set  Xi  ~  iV(0, 1),  Ri  ~  exp(l),  and  n  =  3000.  For  investigation  times  9  = 
.2,  .4,  -6,  .8,  and  1,  we  construct  an  approximate  95%  confidence  interval  for  giNviu*j^y) . 
The  results  are  summarized  in  Table  1.  As  expected,  detection  rate  decreases  as  investigation 
time  increases.  Compared  to  our  result  from  Example  5,  gLoixfo)  =  .295,  the  simulation 
results  indicate  that  reasonably  small  investigation  times  cause  no  significant  degradation 
in  detection  rate.  The  optimal  cycle  length  Un  decreases  with  the  investigation  time — this 
happens  because  investigation  time  causes  more  targets  to  renege  before  detection,  and  the 
perimeter  traversed  is  chopped  at  both  endpoints  to  counter  this  fact.  The  row  headed  by  a„ 
shows  the  optimal  origin,  as  determined  by  the  simulation  and  optimization  results;  a„  does 
not  exactly  equal  —Un/2  due  the  finite  number  of  regenerative  cycles  generated  during  the 
simulation. 


6  Final  Remarks 

In  this  paper  we  developed  models  of  border  surveillance  when  inhltrators  arrive  according 
to  a  Poisson  process  and  can  renege  after  a  random  amount  of  time.  For  both  leap-to-origin 
and  back-and-forth  trajectories,  we  derived  an  expression  for  the  detection  rate  given  the 
reneging,  arrival  location  distribution  and  the  length  of  the  patrolled  sector  (the  decision 
variable).  We  extended  this  analysis  to  study  the  effects  of  sensor  velocity  and  investigation 
time. 

Our  analysis  suggests  that  there  is  not  a  signihcant  difference  in  the  target  detection 
rate  among  the  two  trajectories  considered  when  the  cycle  length  is  very  small  or  very  large. 
When  each  target  is  investigated  upon  detection,  we  hnd  that  the  cycle  length  that  maximizes 
detection  rate  decreases  as  investigation  time  increases,  in  order  to  compensate  for  the  extra 
reneging  of  targets  that  occurs. 
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A  basic  limitation  of  our  models  is  that  the  sensor  footprint  is  one-dimensional.  This 
assumption  is  a  reasonable  representation  of  reality  when  the  size  of  the  sector  is  significantly 
greater  than  the  sensor’s  footprint  width. 

Given  that  our  modeling  assumptions  are  relatively  mild  and  the  operational  scenarios 
we  describe  are  currently  being  employed  (as  in,  for  example,  the  U.S. — Mexico  border; 
see  [1]  and  [17]),  we  believe  that  the  sequence  of  sensor/target  conhgurations  we  pose  in 
this  paper  yields  valuable  insights  into  the  operation  of  sensor  platforms  in  realistic  border 
control  settings. 
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